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Abstract. 

Within the framework of the mode coupling theory (MCT) of structural relaxation, 
mechanisms and properties of non-ergodicity transitions in rather dilute suspensions of 
colloidal particles characterized by strong short-ranged attractions are studied. Results 
building on the virial expansion for particles with hard cores and interacting via an 
attractive square well potential are presented, and their relevance to colloidal gelation 
is discussed. 
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1. Introduction 

Colloidal gelation is a non-equlibirium transition observed in dispersions where short 
ranged attractions exceeding a few ksT are present [[[], |], |3|. It is accompanied by the 
formation of (denser) domains, which coarsen and finally freeze when the gel line is 
traversed [|], ||. The resulting gels are amorphous solids exhibiting finite elastic moduli 
0, 0]. Often, the gel structure is fractal [Q, [|, typically at low densities ||, and often low 
angle scattering peaks are observed, which are reminiscent of spinodal decomposition; 
in this case gelation manifests itself by an arrest of the time dependence of the peak 
position [0, ||, [| [TT], M. This explicit dependence on the time since quenching the 
suspension shows that colloidal gelation is a non-equilibrium phenomenon [[J . 

Colloidal gelation has been observed in colloid polymer mixtures j|, [|] , in solutions 
of sterically stabilized colloidal particles if the solvent quality is decreased || [| |7J , in 
emulsions H , in solutions of charge stabilized particles upon changes of the salt content 



and in protein solutions TTTJ] . In the last case, it prevents protein crystallization 
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and thus is an important obstacle for the collection of protein structure information 



H PL P 



Previous theoretical explanations of colloidal gelation have focussed on the 
possibility of a dynamic percolation within a gas-liquid phase coexistence region which 
is metastable with respect to gas-solid coexistence [[| . The gel is assumed to form when 
the largest cluster spans the sample and the small angle scattering peak consequently 
arrests. Long ranged density fluctuations could arise from quenching below metastable 
(Hl5|) spinodal lines as argued in Refs. 0, || ||, [lOfl . Fine-tuning of the attraction 



strength allows for the observation of all aspects of spinodal decomposition in this 
systems without arrest and gelation || §]. At low densities irreversible cluster-cluster 
aggregation provides a mechanism explaining the observed small wave vector structures 



17| . A review of experiments and theoretical approaches up to 1997 has been given by 
Poon and Haw [|J. 

Recently an alternative explanation for the arrest of the dynamics at colloidal 



gel transitions has been put forward |L8L 19|. Non-ergodicity transitions triggered 



by the local dynamics were suggested. Solutions of equations of the mode coupling 



theory (MCT) for structural relaxation ^TJ were studied which exhibit a localization 



transition explaining a number of properties of colloidal gelation. At low densities but for 
strong short-ranged attractions, particles are tightly bound to (asymptotically infinite) 
clusters and large elastic moduli result |L8], |19[. In this limit, the MCT-equations can be 
simplified to one-parameter models which capture a number of the pertinent physical 
mechanisms although this neglects other aspects of the obtained colloidal gels which are 
connected to their self-similar structure on mesoscopic length scales. 

In the following, the asymptotic models connected to colloidal gelation shall be 
discussed further, where, especially, the existence of a divergent cluster size and its 
related dynamics is of interest. The strength of the required attractions at gelation is 
estimated also. Additionally, two unphysical aspects of the previously obtained MCT 
results for attractive Yukawa potentials and for Baxter's adhesive hard sphere model 



(AHS) |g2[ [18| are reconsidered using the controlled low density virial expansion [23 



Adopting a very common model in colloid science, hard spheres interacting with a square 



well potential are studied p4 
2. Basic equations 

The most simple MCT equations for the normalized time-dependent intermediate 
scattering functions at wave vector q, & q (t), shall be studied, as appropriate for the 
description of the structural relaxation of colloidal suspensions (characterized by a 



Brownian short-time diffusion coefficient D s q ) |25|, 20|, |2l 



d t <S> q (t) + q 2 D s q {<S> q (t) + f dt'm q (t - t') d t ,%(t')} = . (1) 

J 

Here, the generalized longitudinal viscosity is given by a mode coupling functional, 
m q {t) = Fqdc), [$(£)]), which is uniquely specified by the equilibrium static structure as 
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given by the direct correlation function c q : 

^([c], [/]) = ^ / ^ 5(k+p = q) (q • pc p + q • kc fe ) 2 S fc S p / fc / p . (2) 
where p is the particle number density and c q determines the static structure factor via, 

Sg = 1/(1 - pCq). 

Non-ergodicity transitions of MCT dynamical equations are obtained as bifurcation 
points of the non-linear equations for the g-dependent non-ergodicity or Edwards- 
Anderson parameters f q , which, in the idealized MCT, are defined as $ 9 (£ — > oo) = f q : 

T^ L T = FM,[f})- (3) 

Whereas vanishing non-ergodicity parameters, f q — 0, indicate fluid or ergodic states, 
finite ones, f q > f q > 0, spring into existence at larger interactions and signal the 
arrest of density fluctuations into non-equlibrium states, which are amorphous solids 
like gels or glasses. The idealized MCT assumes a separation of time scales describing 
the arrested structures as truly non-ergodic states, although further slow transport 
processes may lead to a decay of the arrested structures at much longer times; see the 
reviews [2C. 21] for discusions and extended theoretical approaches. 

Whereas glass transitions have been studied extensively and are present in the 
studied colloidal suspensions at higher packing fractions also [f|, non-ergodicity 
transitions into gel states exhibit novel phenomena and anomalies connected to the low 
packing fractions and the resulting tenuous ramified gel structures. The experimentally 
observed small-angle scattering peaks indicating domain formation after quenches into 
states far from equlibrium are one aspect, which however is neglected here, because short 
local length scales are considered where local equilibrium may be assumed. Another 
aspect, the existence and growth of tightly bonded clusters is signalled by the local 
dynamics; see below. 



3. Virial expansion input 

Specifying the static input, i.e. c q , the bifurcation lines and the resulting long-time 
dynamics can be determined as functions of the thermodynamic control parameters, viz 
packing fraction, 0, and temperature, T. For colloidal particles characterized by steric 
repulsion and short-ranged attraction, the square well potential is widely used PI] : 

!oo < r < a 

-U a < r < cr + A (4) 
a + A < r . 

The virial expansion provides a controlled approximation for dilute systems and 
determines the direct correlation function for low packing fractions as: c(r) = f(r) + 
0((p), where / = e~ u<yr ^ kBT — 1 is the Mayr cluster function [p3]| . The gel transitions 



can be studied in an asymptotic model of vanishing packing fraction and temperature, 



Colloidal gelation and non-ergodicity transitions 



4 



which describes a dilute system of strongly interacting particles |T8| , |I9 |. Within the 
virial expansion this limit is obtained for: 

12A<f)A 2 

(p — > and A — > oo, so that T v = - = constant , (5) 

IT (T 

where A = e u °/ kBT —l. Note that this results in S q = 1, because of <fr — > 0, indicating that 
structural correlations can be neglected in explaining the gel non-ergodicity transitions. 
Thus, in this limit, glassy dynamics, due to caging at higher packing fractions, is absent 
as well as dynamics connected to critical phenomena, where long-ranged correlations 
are important; note that an increase in the compressibility (or S q= o) would require lower 
temperatures than considered in Eq. ([5]) as the virial estimate indicates: S q ^o oo for 
(j)AA/a — 1/24 + 0((j), — ). Moreover, as expected on physical grounds, only attractive 
potentials can lead to the limit Eq. and thus to gels, because A would remain 
bounded for pure repulsions; this is violated by the mean spherical approximation used 
in 



18, 19 



4. Results for the gel structures 

The results for the gel form factors f q are of immediate interest as they provide 
information on the spatial correlations of the arrested solid-like structures. 




Figure 1. Critical non-ergodicity parameter of strongly interacting dilute fluids 
described by Eq. (||) for the well sizes A/er = |, |, 1 and in the limit A — > 0; most 
curves are not shown below qA = 0.3 because of numerical inaccuracies. The inset 
shows the variation of the critical coupling parameter T£ with well size A. 

In the specified limit of a strongly interacting dilute colloidal fluid, the mode 
coupling functional simplifies to a two parameter model as its density and temperature 
dependence enters via T v from Eq. (|j) only: .F ? ([c], [/]) — > ^(T^, [/]). Figure 1 
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shows non-ergodicity parameters at the transition points for three different well sizes, 
A/ a = 1, |, |, and in the limit A — > 0. At the gel transition, bonds of a length given 
by the attraction range are formed between the particles; the well size A sets the length 
scale for In the non-ergodic states, fast particle rearrangements are possible on 
shorter distances, causing the decay of <& g (£ = 0) = 1 down to f q , but small wave vector 
fluctuations are progressively suppressed upon decreasing A. Only for bond-lengths or 
well sizes A > 0.1, appreciable large distance temporal fluctuations of the gel structure 
exist. 

The additional limit A/ a <C 1 is of interest, because the asymptotic model defined 
by Eq. (|]) further simplifies and the results from Fig. 1 indicate that this limit 
qualitatively captures the gel structures (i.e. ff) for A below ca. 0.1. Moreover, the 
elastic constants of the gel increase strongly in this limit of short-ranged attractions, 
and thus observation of the non-ergodic gel states is more likely, as will be argued below. 




Figure 2. Non-ergodicity parameter (full red line) and critical amplitude h q (long 
dashed red line) at the critical point of the virial model for the square well potential. 
The thin short dashed black, green and blue lines show the corresponding results from 
the AHS model using the wave vector cut-offs: k max a = 40, 80 and 160; the effective 
width is found as A^ s = 3.68/fc max - 

A naive application of A — > leads to the AHS-virial result: c q — > An A Acr sm g g<J . 
Because of its slow algebraic decay for q — > oo, the wave vector integration in Eq. (|2|) 
does not converge and the results depend on the chosen wave vector cut-off k mSuX . In 



the AHS-model this holds for all temperatures and packing fractions as argued in fl8 
However, first entering the direct correlation function into Eq. and then performing 
the limit A — > leads to a different and almost everywhere convergent memory kernel: 



Wv, if}) = i£ d k (^) (i - cos *0 h /la-ki (6) 
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where q = qA denotes the rescaled wave vector. Tight localization of the particles 
is predicted by the solutions to Eqs. ([$]) and (§) shown in Fig. 2, because the q- 
width of the Edwards-Anderson parameters f q is given by the inverse of the narrow 
well width A. Binding to in this limit infinite clusters is described because the wave 
vector dependent longitudinal modulus diverges like l/q 2 , even though it is connected 
to the total force acting among all particles, q 2 T q oc (F* ot (t — > oo)F tot .(0)) for q — > 0, 
where the total force among all particles vanishes according to Newton's actio-reactio 
principle, F tot . = 0. In the specified limit, however, the particles experience forces from 
infinitely removed particles belonging to the same cluster, and thus q 2 T q stays finite as 
the true F to t cannot be determined. Collective and single particle density fluctuations 
consequently become identical, $ 3 (t) = & q (t), leading e.g. to the limit f q — > 1 — (qf s ) 2 
for q — > 0, where the particle root mean square displacement, r s , approaches (half) the 
well width, f c s = r c s / A = 0.48, at the transition point; see also |L8|, [L9| . 

Figure 2 shows the solution of Eqs. (0) and @ (evaluated with 400 grid points 
spaced at 5q = 0.025) at the critical point = 1.42 . . ., which corresponds to rather 
modest attractions; e.g. for A/a = 1/100 and (p = 0.1 one finds U /k B T = 3.56. 
Solutions from the AHS-virial approximation are also included for three different cut- 
offs k max o~ = 40, 80 and 160. One notices that for small wave vectors the solutions can 
be scaled onto the square well result; also the AHS critical coupling parameter scales as 
one would expect: T£ AHS cx l/k metK . However, as all wave vectors for example enter the 



determination of the exponent parameter A [^6| [2(], appreciable differences result: 
A^hs = o.65 compared to A^ w = 0.79. 

5. Results for the dynamics 

The bifurcations to non-ergodic solid-like states in Eqs. (|3|) and @ also lead to 
anomalies in the long-time dynamics which provide the most detailed information about 
the gel-formation and the connected transport mechanisms. We focus on the limit 
of narrow attractions, A <^ a, as clustering is specific to the considered low-density 
transitions, and it is thus of interest to study the consequences of the small-g divergence 
of the longitudinal modulus. 

In the limit given by Eq. ([5]) and for A < u, introduction of a rescaled time, 
t = tDo/A 2 with Dq = D s ao , eliminates the transient parameters and leads to 
simplified equations of motion for the coherent (and identically for the incoherent) 
density fluctuations: 

d~ t <S> q (t) + q 2 § q {) + T dt' mq(t - t') d v $- q (t') = , (7) 
Jo 

where the friction function is given by: m q (t) = g 2 jF^(r^, [/]), and the initial variation 
is $ ? (t) = l-q 2 t + .... 

Figure 3 shows intermediate scattering functions for a range of interaction 
parameters and wave vectors. The two-step relaxation can be analyzed with the 



techniques developed for MCT-glass transitions |20|, It exhibits the factorization 
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Figure 3. Intermediate scattering functions &q(t) determined from (f7|). For the 
interaction parameters T v /T% = 0.9, 0.99, 0.999, and 1.01, corresponding to the short- 
, long-dashed, solid and dashed-dotted curves, solutions for the wave vectors q = 0.8 
(red), 1.6 (blue) and 3.2 (green) are shown. The thin horizontal lines indicate the 
critical non-ergodicity parameters f~ which give the long-time limits, f~ = &q(t — > oo) 
at T£ = 1.3647 for the chosen discretization, Sq — 0.08 and 100 grid points. 



property during an intermediate time window: $q(t) — f~ + hqG x (t) for |$q — f~\ <C 1, 
which enlarges upon approaching T£. The /3-correlator G x , captures the sensitive 
dependences on time and interaction strength, and the amplitudes /- and hg (see Fig. 
2) describe the gel structure and the spatial correlations of the mechanism dominating 
the bond-formation. Power-law decays around the plateau value f c and the divergence 
of a first scaling time are contained in G x (t). 

The near arrest of the <&q(t) around /| leads to an intermediate elastic behavior 
of the gel, which in the studied limit follows from the asymptotic result for the shear 
modulus G: 

7j-F /t- 2 roo 

G(t) = (k B Tp) j o dk (1 - cos k) . (8) 

Figure 4 shows the corresponding storage and loss moduli which, in the fluid, exhibit a 
near elastic plateau with power-law corrections predicted by the /5-correlator G x . Using 
the initial variation of a high-frequency divergence, G'(u) oc G"(u) oc Cj 1 ^ 2 for 

u — > oo results from the integration in Eq. (^|) at high k. This result, which is familiar 
for colloidal hard spheres in the absence of hydrodynamic interactions [23|, cannot be 
seen in Fig. 4 because of the chosen small wave vector cut-off. 

Importantly, Eq. (H) predicts that large, (9((cr/A) 2 ), elastic moduli result upon 
formation of the gels for V > if the particles interact via short-ranged attractions 



27fl . The elastic modulus at the transition, G c , which follows from Eq. (J8|) and the 



critical gel structure factors, /|, is included in Fig. 4. 
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Figure 4. Storage, G'(uj) (blue), and loss, G"(u>) (red), moduli according to Eq. (^J) 
for the interaction parameters and line styles of Fig. 3. The horizontal line indicates 
the critical elastic modulus, G c , of the gel at the transition, r„ = T c v . 



The final decay of the intermediate scattering functions or of the shear modulus G 
is described asymptotically by the second MCT scaling law, which entails the existence 
of a second set of divergent time scales and of r„-independent and non-exponential 
relaxation functions |2"T|| . The small-g divergence of the longitudinal modulus leads 
to a diffusive behavior of the final relaxation times, fg oc T?°(r£ — r„)~ 7 , with the 
I^-independent amplitude f|° oc l/q 2 . Here the exponent 7 is determined by A, and 
follows as 7 = 2.85. The diffusive behavior is apparent in Fig. 3 and stresses that the 
final melting of the precursor gel structure in the fluid phase requires transport over 
large distances which scale like a/ A 2 as follows from Eq. (^|). Again this indicates the 
presence of infinite clusters which arrest and cannot melt anymore at and above the gel 
transition at IT,. 



6. Conclusions 



We have provided further support for the suggestion that colloidal gelation is caused by 
non-ergodicity transitions triggered on local distances |[8|, [L^ . 



It was shown that non-ergodic structures characterized by strong elastic moduli 
are obtained in dilute solutions of colloidal particles interacting with short-ranged 
attractions upon lowering the temperature. Bonds with a length given by the attraction 
range are formed between the particles if the strength of the attraction increases to a few 
kgT. These values are in qualitative agreement with phenomenological considerations 
for reversible gelation |T0L [^]. We suggest, that this formation of long-lived bonds is the 
rate limiting step for (transient) gelation in colloidal suspensions. 
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The formation of (percolated) clusters is indicated by the long-ranged force 
correlations of the arrested structures which cause the small wave vector divergence 
of the mode coupling functional T q . The description of the mesoscopic gel structure 
and of the related domain coarsening kinetics, however, has not yet been incorporated 
into the present approach. 

The long-ranged force correlations at the gel non-ergodicity transitions cause 
that collective and single particle dynamics become identical for long times, and 
that one mode coupling functional determines all time dependent small wave vector 
quantities, like mean square displacement or moduli; such connections have already 
been exploited to measure the elastic response via light scattering techniques |28[ . 
Moreover, a diffusive final relaxation of the colloid density fluctuations results close to 
the gel transitions, r q oc l/q 2 , which suggests that a generalized Gaussian description, 
$ 9 (t) = exp { — (q 2 /Q)Sr 2 (t)}, as suggested by Segre and Pusey p9j , can provide a 
reasonable description of the long time dynamics for all wave vectors; note that for 
colloidal suspensions at the glass transition this ansatz failed at small wave vectors 
because of the non-diffusive structural relaxation [^, [3(J . 

Using the exact low density virial expansion, two at first sight unphysical aspects 
of previous MCT calculations were identified as artefacts of the approximate static 
structural input. First, the cut-off dependence found within the AHS model |22, 
was analyzed for low densities. The interesting conclusion is that MCT predicts the AHS 
model to be in the non-ergodic state for any non-zero particle concentration. Second, 
it was clarified that gelation requires attractive interactions. 

Dynamic light scattering measurements as preformed by van Megen and coworkers 
at the colloidal glass transition [[H], [!2|, [£| or further viscoelastic measurements |7|, ^] 
can provide crucial tests of our approach. First measurements of the non-ergodicity 
parameters support it [34|. Dynamic light scattering measurements of Krall and Weitz 
|35| also show intriguing connections and require further study. 
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